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ABSTRACT 

Quasar accretion disks are believed to form stars by self-gravity. Low Luminosity 
Active Galactic Nuclei (LLAGN) are much dimmer galactic centers, and are often 
believed to be quasars that ran out of gaseous fuel. LLAGN accretion disks should 
thus co-exist with thousands to millions of stars or proto-stars left from the previous 
stronger accretion activity. In principle, these stars may produce several important 
effects: (i) contribute to the optical/UV spectra of some LLAGN; (ii) reprocessing of 
the stellar radiation in the dusty disks could dominate the LLAGN infra-red spectra; 
(iii) deplete the (accretion) gas disk much faster than it can accrete onto the supper- 
massive black hole (SMBH); (iv) stars, individually or in groups, may slow down and 
modulate the accretion flow significantly due to their inertia. In this way they may 
produce the LLAGN cut-off disks; (v) alternatively, frequent enough stellar collisions 
and resulting stellar disruptions could keep the inner disk empty. Here we explore these 
ideas. We find that, despite "low" luminosities of LLAGN, unrealistically high stellar 
densities are required to make a sizable radiative contribution to the (HST) optical/UV 
spectra of these galactic nuclei. Stellar contribution to the infrared spectrum is more 
likely. Further, if LLAGN are in a quasi steady-state for as long as 10 7 years or more, 
too high stellar densities would again be required to significantly affect the dynamics 
of accretion flow. However, if LLAGN are "short" -lived phenomena, e.g. t ,$ 10 5 years, 
quiescent states of quiescence-outburst cycles, then embedded stars may be much more 
important through the mass effects (iii) - (v) . With observations of LLAGN becoming 
progressively better, it will be more and more difficult to neglect the presence of close 
stars in and around nuclear accretion disks. 



1 INTRODUCTION 



LLAGN are galactic centers intermediate in luminosity and 
in some key properties between quasars (bolometric lumi- 
nosity Lb ~ 10 46 erg s _1 ) and inactive galactic centers 
(IGC), such as our own (Sgr A*), whose luminosity is about 
L b ~ 10 36 erg s" 1 . LLAGN are also known as LINERs (low 
ionization nuclear emission-line regions; Heckman 1980) or 
as "dwarf" Seyfert nuclei (e.g. Ho, Filippenko & Sargent 
1997, 2003). Typical nuclear bolometric luminosities of these 
sources are Lt, 10 42 erg s _1 (e.g. Ho 2003). It is well known 
that most present day galaxies harbor supper massive black 
holes (SMBHs), implying that LLAGN and IGC too had to 
be bright accreting sources at some earlier times. LLAGN 
and IGC are thus quasar relics or quiescent phases of a long- 
term evolution of galactic nuclei. Some of the processes op- 
crating in the quasar phase may have long-lasting effects 
and continue to influence SMBH activity in the present dim 
epoch. One such process is star formation inside quasar ac- 
cretion disks due to self-gravity at distances greater than 
~ 10 3 Schwarzschild radii (e.g. Paczyriski 1978; Kolykhalov 
& Sunyaev 1980, Shlosman & Begelman 1989; Goodman 
2003), where Rs = 2GMbh/c 2 , the Schwarzschild radius 
for the black hole with mass Mbh- The goal of this paper 



is to investigate some of the possible effects of the presence 
of these close and fast moving stars on the LLAGN phe- 
nomenon. 

Let us summarize the observational and theoretical 
framework in which bright AGN, LLAGN and inactive 
galactic centers are currently understood. Quasar accre- 
tion disks are believed to be the "standard" accretion disks 
(Shakura & Sunyaev 1973), possibly supplemented by an 
X-ray emitting corona. The cold accretion disk apparently 
extends down to the last stable orbit, as evidenced by 
a very prominent thermal-like UV feature called the Big 
Blue Bump (BBB) and by the broad Fe Kq lines observed 
in many cases (it is our personal opinion that the non- 
detection of broad Fe Ka lines in other sources is caused by 
an over- ionization of the accretion disk surface). The near- 
Eddington luminosities of quasars, short time scale variabil- 
ity, and the high radiative efficiency with which most of the 
SMBHs are inferred to have grown (Yu & Tremaine 2002) 
are also best explained in the context of the standard ac- 
cretion disks since they indeed have very high radiative effi- 
ciency. 

LLAGN, and our Galactic Center (GC) black hole, 
Sgr A*, are distinctly different in this regard, possibly re- 
quiring a radiatively inefficient accretion. Chandra X-ray 
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observations constrain the amount of hot gas at the capture 
radius in Sgr A* and M87 (Baganoff et al. 2003; Di Matteo et 
al. 2003), which allows one to estimate the Bondi accretion 
rate onto the SMBHs. For other sources where direct X-ray 
data are not constraining enough, Ho (2003) notes that the 
centers of the galaxies hosting LLAGN should be filled with 
hot winds from young stars (similar to our GC), yielding 
an accretion rate estimate. In addition, the accretion rate 
through a cold accretion disk may be estimated by using 
the optical/UV data in some LLAGN (e.g. Quataert et al. 
1999; Ptak et al. 2004). All of these methods lead to the 
same conclusion - the expected luminosity from a thin stan- 
dard accretion disk fed at the deduced accretion rate is some 
1 — 4 orders of magnitude higher than actually observed. 

In addition, in stark contrast to quasars and Seyfert 
galaxies, LLAGN do not show the BBB in their spectra 
(Ho 1999), have no short time scale X-ray variability and 
no broad Fe Ka lines (e.g. Dewangan et al. 2004; Ptak et 
al. 2004). This strongly argues for the absence of the cold 
accretion disks in the inner ~ 10 2 — 10 3 i?s. Theoretical in- 
terpretation of this fact is that the cold disk "evaporates" 
at few hundred to few thousand Rs (e.g. Quataert et al. 
1999; Fig. 2 in Ho 2003; Ptak et al. 2004) to form an in- 
ner hot Non-Radiative Accretion Flow (e.g. Narayan & Yi 
1994; Blandford & Begelman 1999). The most frequently 
discussed evaporation mechanism is that of Meyer & Meyer- 
Hofmeister (1994). If the ions of the hot flow are much hotter 
than the electrons, then the flow is radiatively inefficient and 
is much less luminous than the standard flow, explaining the 
discrepancy. In the least luminous (in terms of Lb/^Edd, 
where I/Edd = 1-3 x 10 46 Ms erg s _1 , the Eddington limit 
for the black hole mass Ms = Mbh/10 Mq) cases, such as 
Sgr A* , it is possible that a cold accretion disk does not even 
exists (e.g. Falcke & Melia 1997; Narayan 2002). 

The best studied example of an inactive galactic center 
is the center of our own galaxy, Sgr A*. Genzel et al. (2003; 
see also Levin & Beloborodov 2003) discovered that most 
of hot young (t < 10 7 years) stars in the GC lie in one of 
the two stellar disks, strongly suggesting these stars were 
formed by the self-gravity in two massive accretion disks 
existing in earlier times. In the case of the GC, these young 
stars actually dominate and determine the SMBH activity, 
since their radiative output is some ~ 5 orders of magnitude 
greater than that of the accretion flow, and since the hot 
winds produced by these stars are believed to be the main 
source of fuel for Sgr A* (e.g. Coker & Melia 1997). 

It is thus likely that important interactions between 
close stars and the accretion flow may be taking place in 
LLAGN as well. Star formation may even continue in the 
LLAGN phase as well (e.g., see Fig. 1 in Goodman 2003; 
also Milosavljevic & Loeb 2004). Quite generally, close stars 
should assume more prominent roles in the LLAGN (and 
IGC) phase than they do in quasars because the mass and 
the luminosity output of the gaseous accretion disk are mod- 
est in these sources. The possible effects of the close stars 
can be divided into two categories: the luminosity and the 
mass effects. The former are due to the combined stellar 
luminosities; whereas the latter are due to stellar mass - 
dynamical and accretion effects. We find that the luminos- 
ity effects, besides reprocessing of the stellar light into the 
infrared region, are unlikely to be important due to existing 
observational constraints. The mass effects are much more 



likely to be important especially if the disk inner regions are 
"light" as is the case when the disk is unstable to a limit 
cycle instability (e.g. Siemiginowska, Czerny, & Kostyunin 
1996; Lin & Shields 1986) or is emptied out in the inner part 
for other reasons. 



2 SELF-GRAVITY IN QUASAR DISKS 

We shall not discuss the implications of self-gravity for 
quasar disks because these have been investigated previously 
by a number of workers (e.g. Paczyhski 1978; Kolykhalov & 
Sunyaev 1980; Shlosman & Begelman 1989; Goodman 2003). 
Quasars are believed to be accreting near the Eddington 
limit and hence their accretion disks need to be quite mas- 
sive. The disks are expected to form stars at radii for which 
the Toomre (1964) parameter satisfies 

_ Cs n _ n 2 H _ hmbh „ m bh h 

Q ttGE ~ ttGE tyER 3 M d (R) R ( ' 

where c s is the isothermal sound speed of the disk, O is the 
Keplerian circular rotation frequency, S(i?) is the disk col- 
umn depth at radius R, H is the disk half-thickness (pressure 
height-scale), and finally, Md(R) = T,nR 2 is approximately 
the disk mass within radius R. Following Goodman (2003), 
the radii which satisfy this conditions in the standard accre- 
tion disk are 

/ -2\ 2 /9 

i?>i? Q=1 ^2.7xl0 3 i? s (^gU , (2) 

where a_i is the normalized Shakura-Sunyaev viscosity pa- 
rameter, m — M/M E dd, where M E dd = -^Edd/ec 2 , with ra- 
diative efficiency e = 0.1. In this paper we shall use the 
notations in which a_i = a/10 -1 , rh_3 = m/10 -3 , etc. 

It has been argued that, in addition to Q 1, a second 
condition for star formation by self-gravity in accretion disks 
must be fulfilled: the cooling time should be shorter than a 
fraction of the Keplerian rotation period (e.g. Shlosman & 
Begelman 1989; Gammie 2001). However, recent simulations 
by Johnson & Gammie (2003) demonstrated that the insta- 
bility domain is much larger. In the non-linear instability 
stage, the local cooling time may be orders of magnitude 
smaller from that found in the unperturbed disk model. 
Therefore Johnson & Gammie (2003) concluded that it is 
extremely difficult to prevent even the marginally stable, 
Q ~ 1, optically thick disk from fragmenting and forming 
stars. Since the formation of gravitationally bound struc- 
tures is extremely fast (Toomre 1964), it is then feared 
that most of the disk mass would be used up to build 
up stars rather than feed the SMBH. Sirko & Goodman 
(2003) demonstrated that the amount of heat required to 
stop star formation in the disk outer regions is excessive 
(i.e. larger than Z/Edd for a typical quasar), and moreover 
the observed quasar SEDs rule out such a heating. Thus it 
seems inevitable that stars will form in quasar disks at radii 
Si 10 3 Rs. Having formed there, the stars can migrate (in- 
spiral) radially through exchange of the angular momentum 
with the disk (e.g. Syer, Clarke & Rees 1991), so their radial 
distribution may evolve significantly on time scales that are 
long enough. 
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It is not clear how much mass would be converted to 
stars in a single accretion outburst episode, but if we assume 
that this is a fraction q < 1 of the disk mass, then obviously 

M„N, = qM d w 1 |m bh ~ 3 x 10 5 M Q 1 M s , (3) 
Q R Q 

where we set H/R ~3x 10 -3 (see equation N* and M* 
are the total number of stars created and a typical stellar 
mass, respectively. It is also unclear how exactly this stellar 
distribution evolves with time under the interaction with 
the accretion flow, since stars may be swallowed or tidally 
disrupted by the SMBH if they migrate close enough (Syer 
et al. 1991), they may explode as supernovae (e.g. Syer et 
al. 1991; Rozyczka, Bodenheimer, & Lin 1995), destroyed by 
direct star-star collisions, and evolve under the influence of 
long-range two-body interactions (e.g. Bahcall & Wolf 1976). 

It is outside the scope of this paper to discuss ways for 
quasars to maintain high enough mass transfer rates into the 
SMBHs when the self-gravity sets in their accretion disks. 
Now, suppose that the quasar mass supply has substantially 
decreased, so that most of the gaseous disk mass has been 
used up, and that some number of the "original" self-gravity 
created stars is still caught up in the disk. The quasar pre- 
sumably becomes a LLAGN. What would then be the role 
of such stars? 



T 
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Figure 1. M81 infrared to UV SED compared to the spectrum 
of the central parsec of the Milky Way (solid curve). The dotted 
line shows the spectrum of a nuclear star cluster that "fits" the 
M81 data. Note that the latter model requires too high a spatial 
density of stars and hence it is ruled out on physical grounds (see 
text). 



3 DIRECT CONTRIBUTION TO LLAGN 
SPECTRUM 

Even with the HST, except for very few nearby galaxies, 
galactic centers are not resolvable to scales better than ~ 1 
pc in the optical/UV. Stars formed in a self-gravitating disk 
at a distance of R = 0.04(A/ B h/4 x 10 6 M )r 5 pc, where 
r = R/Rs, and r 5 = r/10 5 , cannot be directly separated out 
from the accretion disk spectrum. In addition, the galactic 
centers are viewed through a large column depth of galac- 
tic stars (i.e. the galactic bulges). However, Filippenko & 
Sargent (1985), and Ho, Filippenko & Sargent (1997), cre- 
ated a large template of galaxies without or with very weak 
emission lines. Using a sophisticated technique to fit the ab- 
sorption lines in the spectra of LLAGN, they were able to 
subtract out contribution of galactic stars, and to obtain the 
nuclear (very weak compared to the galaxy) continuum and 
broad emission lines of LLAGN. 

The stars that we are interested in are very close to 
the SMBH, and their circular Kepler velocities are vk(t) w 
2000 r~ 1/2 « 2000 (0.1 pc M s /R) 1/2 km/sec, where r 4 = 
r/10 4 . Since the galaxy template contains spectra with 
FWHM velocities of < 500 km/sec (Ho et al. 1997), it is 
most likely that the Doppler-broadened spectra of the in- 
ner stars in the LLAGN would be counted as the continuum 
emission (Ho 2004, private communication). Can such stars 
make a non-negligible contribution to the LLAGN nuclear 
spectrum? Two cases exist. 

3.1 Hot massive young stars 

The stellar spectrum of the inner central parsec of our GC 
is dominated by hot young stars. It is thus interesting to 
compare its spectrum with the better of the available nuclear 
LLAGN spectra. M81 is a nearby LLAGN, with the SMBH 



mass estimated at a value very close to that of Sgr A* (e.g. 
Schodel et al. 2003; Ghez et al. 2003), i.e. M BH = 4x 10 6 M . 
M81 nuclear infrared to UV spectrum is shown as thick data 
points and upper limits in Figure^as compiled in table 2 of 
Ho (1999). Ho (1999) argues against a significant intrinsic 
dust extinction in LLAGN. The spectrum of stars in the 
inner 1 pc in our GC is relatively well known (Mezger et 
al. 1996): there is a contribution of ubiquitous low mass 
stars with effective temperature T c g = 4 x 10 3 K, with a 
combined luminosity of 10 6 Lq, and that of ~ 30 hot young 
stars with the total luminosity of 3 x 10 7 Lq and the effective 
temperature of T e g = 3.5 x 10 4 K. This spectrum is shown as 
a solid curve in Figure^and assumes no extinction (of course 
we actually observe the GC through some ~ 30 magnitudes 
of extinction in the visual). 

Clearly, had there been a similar amount of such hot lu- 
minous young stars in M81 as in the GC, these stars would 
easily out-shine the LLAGN in M81 in the UV band. How- 
ever in the GC the vast majority of the hot young stars 
move with velocities of order £ 500 km/sec, so that the 
subtraction technique of Ho et al. (1997) would have proba- 
bly removed them from the nuclear spectrum. In any event, 
stars moving with so "low" velocities are outside the ra- 
dius where most of the accretion flow physics is played out 
( £ W 5 Rs). The results of Ho (1999) thus clearly show that 
there are no or very few hot young luminous O-stars in the 
inner accretion flow region in most of the LLAGN (LINERS) 
as otherwise they would produce a component similar to the 
BBB. 

Since the history of star formation and evolution in 
an accretion disk is extremely model- and parameter- 
dependent, it is possible that such young massive stars never 
formed in the LLAGN accretion disk. It is also possible that 
these stars have already gone through the supernova phase 
(e.g. Syer et al. 1991), or are on the red giant branch, thus 
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shifting their spectrum from the UV into the optical fre- 
quencies (but see the next section). 



3.2 Low mass stars or evolved red giant stars 

The dotted line in Figure Q is a very simple model that 
"fits" the M81 data much better. In this model the "low" 
mass stars have average temperature of 6 x 10 3 K, with the 
total luminosity of 6 x 10 6 Lq , and just one hot luminous 
star with L = 10 6 L© and T cS = 3.5 x 10 4 K. Clearly, a 
more realistic model with a power-law mass function could 
be designed that would easily fit the ifra-red to UV SED 
of M81 with no contribution from the accretion onto the 
SMBH whatsoever. 

However, there are several arguments against this pos- 
sibility. As we argued above, to appear as AGN contin- 
uum, the spectrum of such stars would have to be strongly 
smeared, say with random velocities of order 2000 km/sec 
or more. For M81, the mass of the stars would exceed the 
deduced SMBH mass, thus leading to a direct contradiction 
with observations. Further, the stellar densities in the inner 
part of the stellar cluster would be as large as 10 14 Mqpc -3 if 
the stars are main sequence stars with M ~ Mq . Such high 
stellar densities are extremely unlikely because stellar phys- 
ical collisions and disintegrations would be expected (e.g. 
Murthy, Cohn, & Durisen 1991). If the stars are much more 
luminous super-giants, then the required stellar number den- 
sity is much lower but because of their large radii, we find 
that the collision time is of the same order, i.e. too short. 
The red giant stars are also easy targets for ablation of mass 
due to star-disk collisions (Armitage , Zurek & Davies 1996) . 
Finally, there is a deficit of such stars in the inner cusp of 
our GC (Genzel et al. 2003). 



3.3 Infrared spectrum: reprocessing in the disk? 

The stellar cusp/cluster could also contribute to heating the 
gaseous disk, especially its surface layers. If the stars are uni- 
formly distributed within radius R, then the effective tem- 
perature of the stellar radiation is 63L\{' i (R/pc)~ 1 ^ 2 Kelvin, 
where L41 is the total stellar luminosity in units of 10 41 erg 
sec -1 . We argued that it is not plausible for the stars to dom- 
inate the nuclear optical/UV continuum of LLAGN, hence 
it seems unlikely that they would significantly add to the 
reprocessed disk emission. 

This may not be the case, however, because the accre- 
tion disk is very thin (equation |5J, and only a very small 
fraction of the luminosity generated by accretion onto the 
SMBH may impact the outer disk. In contrast, the accretion 
disk is geometrically thick compared to stars and would re- 
process a large fraction of stellar light even if stars are ex- 
actly in the disk midplane. If stars are on more general orbits 
(see then roughly half of their radiation is reprocessed in 
the disk, possibly making it important in the mid-infrared 
region where a "bare" disk (without surrounding stars) is 
predicted to be dim. In this case the distribution of the ra- 
diation impinging on the disk will be quite patchy. Therefore 
there may be a broad range of temperatures resulting from 
the stellar illumination of the disk which would obviously 
lead to a broad infra-red feature in the SED of the LLAGN. 
This feature will be spectroscopically much broader than 



that resulting from the disk illumination by the single cen- 
tral source (the accreting SMBH), and it may be also vari- 
able on time scales of weeks to years because of bright stars 
approaching and passing through the disk (e.g. Nayakshin 
& Sunyaev 2003; Cuadra et al. 2003). Stars also heat the 
disk via bodily impacts (e.g. Zentsova 2003, Norman & Silk 
1983) and can increase the effective disk viscosity (Ostriker 
1983). 



4 EMBEDDED STARS AND DISK 
ACCRETION 

Stars orbiting the SMBH in the inner parsec of the accre- 
tion flow will also modify the accretion flow of gas. We as- 
sume that the combined mass of stars in this region is much 
smaller than Mbh, and consider the case when the stellar 
and disk planes coincide exactly. In this limit the relative 
velocity of a star and the gas in the disk surrounding the 
star is negligible. Thus star-disk interactions are qualita- 
tively different. The interaction is best understood based on 
the analogy between this case and that of a better stud- 
ied circumstellar disk with embedded (proto-) planets. If 
the star does not open up a radial gap in the disk (see be- 
low) , then it grows via Bondi-like gas accretion. If the gap is 
opened, the accretion onto the star slows down significantly. 
The star however continues to interact with the disk due to 
gravity (essentially via dynamical friction). In both gap/no- 
gap cases the star can migrate, that is spiral-in, radially. 
Clearly the flows of gas and stars are not independent. 



4.1 Gap formation conditions 

A star co-rotating with the disk interacts with the gas in the 
disk via gravitational torques (e.g. Goldreich & Tremaine 
1980; Lin & Papaloizou 1986a, b). The linear velocity of gas 
at radii greater than the radius at which the star circles the 
SMBH, R, is increased as a result of the interaction with 
the star. Thus its angular momentum increases and it moves 
outward from the star's location. The opposite happens to 
the gas interior to the star's location. Hence the star pushes 
the gas away from its radial position. If this effect is strong 
enough, a radial gap is opened in the disk. 
The Hill's radius R H , 



Ri 



Ah 



3Mbh 



1/3 



R . 



(4) 



defines the star's sphere of influence in the disk: dynamics of 
gas within the sphere is dominated by the star rather than 
the SMBH. 

There are two conditions for formation of the radial gap 
in the disk. The first is that Rr is greater than the pressure 
scale height of the disk, H. Otherwise the gas would be able 
to overflow the star at vertical heights z > Rh- The disk 
height-scale H for gas-pressure standard accretion disk at 
r = R/Rs > 1 is 



H/R = 2.5 x 10~ 3 (aM 8 ) 



-1/10 1/20-1/5 



(•») 



(see Svensson & Zdziarski 1994) . We assume that the accre- 
tion disk is not strongly self-gravitating during the LLAGN 
epoch, i.e. that Q ^ 1. Since 



Close stars in LLAGN 5 



— = 0.26 (aMg) 
Rh 



-1/10^/20^1/5 (MS 



1/3 

— ] <1, 



(6) 



where m* = M,/ Mq. The first condition for a gap forma- 
tion is usually satisfied in LLAGN disks. 

The second criterion is that the rate with which the gas 
is pushed away from the star's location due to gravitational 
interaction is higher than the rate with which the gap could 
be filled in due to viscous gas transport. This is possible 
only when the disk viscosity is low enough. Parameterized 
in terms of the Shakura-Sunyaev a-parameter, this condition 
becomes 



1 / M* \ 2 f R\ 



40 VMbiT; KhJ (7) 

(e.g. Syer et al. 1991). Using the expression for H/R, we get 



4 

3.8 x 10~ 5 ^t-2 



V2„ 



(8) 



It is unlikely that a is as low as this, and hence the low mass 
stars are unable to open a gap in the disk of a "large" SMBH, 
i.e. Ms ^ 1. However, for M B h = 4 x 10 6 Mq, the coefficient 
in equation|H|is a m m = 0.6: low mass stars may open a gap in 
accretion disks of lighter galactic nuclei. Finally, for Ms = 1 
and m* = 10, i.e. a high mass star with M* = 10 Mq, we 
have a m in = 0.38, thus a gap formation is again likely. 



4.2 No gap: accretion onto embedded stars. 
Forming an inner hole? 

Consider now the cases when the stars are unable to open 
up a gap. As shown above, a gap is unlikely to form for high 
mass black holes, i.e. Ms ^ 1, or when the stars are low mass 
stars. Also, gap formation is more problematic at large radii 
i-2 > 1. The stars then grow by accretion of gas from the 
disk (Syer et al. 1991). 



4-2.1 Low mass stars 

For low mass stars, stellar radiation and winds are relatively 
weak. The accretion of gas onto a low mass star is then 
similar to the growth of terrestrial planets in a planetesi- 
mal disk (Lissauer 1987; Bate et al. 2003; Tan & Goodman 
2004). Since Rh > H , accretion onto the star is quasi two- 
dimensional. The accretion rate is determined by the rate at 
which differential rotation brings the matter into the Hill's 
sphere, 



M* = M H ~ AiyRhH pv H ~ 4nR H pCs 



(9) 



where p = E/2/f is the mean disk density. We used the 
fact that the characteristic gas velocity (relative to the 
star) at the Hill's distance from the star, vh, is vh = 
Rn\dQ. J 'din R\ ~ c 3 (Rh / H) since the angular velocity for 
Keplerian rotation is Q. = c 3 /H. Equation^ is valid as long 
as Rh > H since in the opposite case the gas thermal veloc- 
ity becomes important and the accretion would proceed at 
the Bondi rate (e.g., Syer et al 1991). Of course M* cannot 
exceed the Eddington limit value, e.g., 



,Edd 



10 -3 — M Q year -1 
r 



(10) 



where r* =7?*/ Rq, star's radius in Solar radii units. 



4-2.2 High mass stars in the disk 

The enormous radiation field of luminous high mass stars 
pre-heats the disk material at distances ~ Rh from the star. 
If the temperature established is much higher than the un- 
perturbed disk temperature, the flow of gas will be strongly 
affected. Indeed, T eS (R H ) = [L, /AtvgR h ] 1/4 , 



T cS (Rh) = 730 K 



Ll /4 M, 



-1/3 -1/2 



(M t /3OM ) 1 /6 ' 

(where L5 = L»/10 5 Lq) can be larger than the disk effect 
temperature, T c g^: 



(11) 



T cff . d = 110 K 



m_4 



{m 



1/4 



(12) 



The accretion disk next to a high mass star is thus bloated. 
At distances from the star ^ Rh, the star's gravity takes 
over and a thin small scale circum-stellar disk develops. Per- 
pendicular to the disk, strong stellar winds sweep up the 
colder matter. A sort of torus thus forms around the star. 
There will be a thin ionized (Stromgren) layer on the side 
of the torus facing the star. At low disk surface densities 
appropriate to LLAGN, high mass stars loose rather than 
gain mass from the disk; the winds may cool and condense 
on the disk surface 1 . The accretion rate onto the massive 
star is strongly reduced due to the disk pre-heating. The 
disk accretion on the high mass stars is completely negli- 
gible compared with that on the much more numerous low 
mass stars that do not suffer from pre-heating and stellar 
winds. 

The observational appearance of the star will be 
strongly angle dependent. Viewed pole-on to the disk, the 
intrinsic stellar emission is clearly visible, as is the line emis- 
sion resulting from the disk illumination. However, when 
viewed at large inclination angles, the stellar emission is 
hidden from the observer similar to the central regions of 
type 2 AGN. The optical line emission of the inner part of 
the "torus" around the star may still be visible. 

Note that for quasar disks, rh ~ 1 and then 
7cff,d ~ Teff(-Rir) even for the brightest O stars. In addi- 
tion, for smaller radii, i.e. R ~ 100i?s, the disk pre-heating 
is unimportant even for LLAGN. 

4-2.3 Eating the disk from within 

As shown previously, the embedded stars gain mass from the 
accretion disk. Here we consider only the low mass stars that 
are certain to dominate the disk-stars mass exchange due to 
their large number. From the standard accretion theory, the 
gas accretion rate onto the SMBH is 



M B h = A-kRHp\vr\ ~ AirH pac s 



(13) 



since \vr\ ~ a(H/R)c 3 (e.g. Frank, King & Raine 1992). 
Therefore the rate of accretion on the single star compares 
to the accretion rate onto the black hole as 



M H 



R H 



Mbh aH 2 



1 !! 



(14) 



Note that neither the hot stellar wind nor the gas evaporated 
from the disk escape to infinity because the escape temperature 



^ vir 



■ 3 x 10 s Kr7 
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Figure 2. Modified cold disk fits to M81 infrared to UV SED. 
The standard disk with the inner cutoff radius of IOO-R5 is shown 
with the thick solid curve (Quataert et al. 1999 model). The rest 
of the curves assume cold inner disk persisting to the last stable 
orbit but the accretion rate decreasing with radius as prescribed 
(see text for detail). 

Surprisingly, just one star in a LLAGN disk may consume 
the disk quicker than the SMBH, leaving little for the SMBH 
growth. Of course the accretion onto the proto-star is lim- 
ited to M*,Edd> for one. Pre-heating of the gas by the ra- 
diation from the star and the circum-stellar accretion disk 
may decrease the accretion rate further. Nevertheless, even 
if accretion on a single star is very inefficient compared with 
equation [5] a large number of proto-stars in the disk can 
make up for this, and thus the combined mass accretion on 
these stars could well exceed that onto the SMBH during 
the LLAGN phase. 

4.3 Carving out an inner hole? 

If stars embedded inside the disk eat it away before the gas 
can reach the SMBH, then the disk inner region may in prin- 
ciple be almost empty and hence underluminous compared 
with the steady-state standard accretion disk. This could 
then explain the observed lack of the BBB in the LLAGN 
spectra, interpreted as an inner hole in the accretion disk 
(e.g. Quataert et al. 1999). 

We test this idea by comparing the expected disk spec- 
tra to the observations of M81 as an example. Given the 
quality of the data, we model the SED of cold disks con- 
sumed by the stars with a very simple model. We use the 
standard multi-temperature black-body spectrum with the 
exception that the accretion rate, M(R), is a function of ra- 
dius. We neglect stellar irradiation compared with internal 
disk viscous heating. The local disk effective temperature is 
then given by equation 1121 

We assume that stars are only present at radii larger 
than a "self-gravity" radius, R ss . The latter may be smaller 
than Rq=i because of the possible radial migration of stars. 
For simplicity, we choose m(R) — rh(R out )(R/ Rout) p , for 
R sg < R < R ou t, and rh{R) = rn(R sg ) for R < R ag . Figure[2] 
shows the same data as Figure^with several models for the 



disk emission. The solid line shows the alternative cut-off 
disk model similar to Quataert et al. (1999) - the standard 
accretion disk with the inner radius of Ri = 100i?s and the 
disk accretion rate rn(R out ) — 4 x 10~ 3 . The cold flow is 
presumed to evaporate and not to emit in the optical/UV 
region for R < Ri. The dotted line shows the model with 
p = 1, R sg = WORs, Rout = W 5 Rs and m = 4 x 10" 3 . 
Since for p — 1 the energy emitted per a decade in ra- 
dius (oc GMbhM(R)/R 3 x nR 2 ) becomes independent of 
radius, there is a flat part in the infrared frequency part 
(y ~ 10 13 - 10 14 Hz) of the spectrum for this model. The 
model is too dim in the optical/UV. The dashed line shows 
a similar model but with _R sg = 10 3 Rs- The model is clearly 
incompatible with the data. 

Finally, the dot-dashed curve shows a more complex 
model with m = const = 3 x 10~ 2 for 500J?s = Rsg < 
R < Rout, m = m(J? sg )(7?/ii sg ) 6 for 100R S < R < i? sg , 
and rh(R) — rh(R Bg ) for R < R Bg . R ou t = W 5 Rs in the 
model. This model roughly "fits" the data. However, is is 
clearly fine-tuned to have a similar inner radius of WORs 
as the model of Quataert et al. (1999). The "eaten- from- 
within" disk model may fit the data only if the accretion 
rate at the region R ~ few hundred Rs is similar to the cut- 
off disk model, and then drops very sharply. (What hap- 
pens at larger radii is not that important as long as the 
outer disk is not too luminous in the infra-red region). A 
sudden disk depletion at a small range of radii is not im- 
possible if there is a large concentration of low mass stars 
there, but is seems a rather unlikely case to us. The total 
mass of the stars in this disk depletion region would not 
be very large, 8M, ~ M(R out )t ~ 10 3 Af (t/lO 7 years) for 
rh = 4 x 10 -3 . However, the radial size of the 100i?s re- 
gion is only ~ 10 14 cm = 3 x 10~ 5 pc, for the M81 SMBH 
mass, yielding stellar density of order 3 x 10 16 Mq /pc 3 . Any 
star present at such small radii from a pre-existent isotropic 
cluster would have been destroyed in only ~ 3 years. The 
resulting disintegration of the stars would have resulted in a 
very energetic accretion event. Thus this picture would man- 
date periods of quiescence when stars steal most of the mass 
from the accreting flow, and shorter periods of strong accre- 
tion activity onto the SMBH when the stellar collisions or 
supernova explosions (if stars grow to be too massive) release 
the stored gas. As time between these luminosity variations 
is long compared to history of Astronomy, it is not immedi- 
ately clear whether such model is supported by observations 
or not. 

4.4 Gap accretion. Clogging accretion flow with 
embedded stars? 

For higher stellar mass to the SMBH mass ratios, gap for- 
mation occurs (see eqs. |7|&|HJ. The accretion rate of the 
star is reduced by the gap presence, but does not stop com- 
pletely (e.g. Bate et al. 2003) . When the mass of the star is 
small compared to the surrounding disk, the radial motion 
of the star is similar to that of the fluid at the same position 
in the disk. The star thus migrates closer to the SMBH on 
the disk viscous time scale. Because of its small mass, the 
star does not significantly alters the radial surface density 
profile of the disk (Lin & Papaloizou 1986b). While the de- 
tailed disk spectrum may be affected by the presence of the 
star, the effect is quite minor as long as the star's mass is 
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small compared to that of the disk, and the disk is optically 
thick. 

However, when the mass of the star (or any massive 
satellite such as a secondary black hole) is considerable, the 
accretion disk is bound to "notice" the satellite (Syer & 
Clarke 1995). Since the angular momentum exchange con- 
tinues through the gap, the disk interior to the star is able 
to continue accretion onto the SMBH. At the same time, the 
gas at radii larger than the position of the star is blocked by 
the gap. Because of the high inertia of the star, the radial 
velocity of the gas just outside of the gap is smaller than 
its unperturbed (i.e. star-free) value, and the disk surface 
density builds up there (Lin & Papaloizou 1986b; Syer et al. 
1991; Ivanov, Papaloizou, & Polnarev 1999). Syer & Clarke 
(1995) show that the radial migration velocity of the satellite 
becomes 



where vr^ = 3a(H/ R) 2 vk{R) is the disk unperturbed ra- 
dial viscous velocity. If M* S> 7rE_R 2 , when the disk ac- 
cretion is basically temporarily stopped, until the mass be- 
hind the satellite becomes large enough to push it through 
into the black hole. Meanwhile, the inner disk, cutoff from 
the mass supply, empties out and may become quite dim 
(Syer & Clarke 1995). This situation will persist during time 
* ~ M*/M ~ 10 3 m. Af 8 -1 m_ 3 years. This is actually an up- 
per limit, because the overflow of the gap may occur at lower 
values of E, i.e. at somewhat earlier times. After enough 
mass has been accumulated, the accretion resumes, result- 
ing in a bright outburst. The star and the built-up disk 
accrete onto the SMBH on the disk viscous time scale, the 
one that reflects the amount of the stored mass. If the ac- 
cretion rate in such an outburst is high enough, the inner 
disk becomes radiation-pressure dominated and the viscous 
time may be significantly shorter than that of the quiescent 
disk at the same location (this is due to a change in the 
factor (H/R) 2 ). The activity phase is thus relatively short 
compared with the mass accumulation phase. 

Are such blocked accretion flows relevant to the ap- 
parently cutoff disks of LLAGN? If this were the case, the 
LLAGN activity would be periodic, with the LLAGN be- 
ing the quiescent state and a more brighter phase, such 
as Seyfert galaxies nuclei, would be the outburst state. 
The mass of the unperturbed accretion disk, estimated as 
Ma ~ nR 2 'E(R) in the standard accretion flow is 

M d = 60 M s a~ 4/s M 8 u/5 ^ /5 m% 5 . (16) 

The disk is too heavy for a single star to induce serious ac- 
cretion rate modulation for these parameters. For a smaller 
SMBH mass, say 4 x 10 6 Mq, the coefficient in equation ll6l is 
about 0.05 Mq, making a single star a big drag onto the ac- 
cretion flow. The accretion luminosity can then be substan- 
tially modulated. Syer et al. (1991), considering the case of 
a 10 6 Mq SMBH, already mentioned such a possibility. The 
quiescence time can be as long as 10 years in this case. 

Therefore, the possibility of a significant time variabil- 
ity due to presence of stars in the disk, and the "empty" 
observational appearance of the inner disk, cannot be ruled 
out for lower mass black holes. However, the absence of the 
UV bump (BBB) in the LLAGN spectrum is quite a general 
feature, seemingly independent of the SMBH mass. Indeed, 



many of the objects in the Ho (1999) sample have SMBH 
masses higher than 10 8 Mq . The recent observations of Ptak 
et al. (2004) of NGC 3998, a LLAGN with the SMBH esti- 
mated at ~ 10 9 Mq, also shows no clear BBB. Many stars 
embedded in the accretion flow and all drifting together 
could slow down the accretion in the inner regions of high 
SMBH accretion disks, but one would naively expect the 
stars to be smoothly distributed in the disk. Hence embed- 
ded stars cannot explain cutoff disks in high mass LLAGN. 
Therefore clogging up the disk flow with embedded stars 
is possible, but it cannot be the reason behind the missing 
inner accretion disk of LLAGN. 



5 STARS AT HIGH INCLINATION ANGLES 
TO THE DISK 

So far we have been concerned with stars embedded in the 
disk. In general, however, there will be stars on highly in- 
clined (to the disk) orbits, brought in by relaxation of the old 
central star cluster (e.g. Bahcall & Wolf 1976), or from par- 
tial capture of stars by the accretion disk (Syer et al. 1991), 
etc. If stellar orbital planes and the disk plane are signifi- 
cantly mis-aligned, then most of the interactions take place 
when the stars hit the disk and pass through it (e.g. Ostriker 
1983; Norman & Silk 1983; Zentsova 1983; Perry & Williams 
1993; Nayakshin et al. 2004). The stars move through the 
disk at very high Mach numbers and produce strong shocks. 
Quite generally, the stars are expected to enhance (speed- 
up) the accretion flow because of the additional angular mo- 
mentum transfer mechanism (Ostriker 1983) and the extra 
heating. Thus this effect works in the opposite direction to 
that discussed in H4.4I 

For radii smaller than some value Ri, the angular mo- 
mentum transfer in the disk is dominated by star-disk col- 
lisions rather than by the internal disk viscosity. The value 
of -Ri is a function of the number density of stars in the in- 
ner part of the stellar cusp and their orbits (see Nayakshin 
et al. 2004 for estimates for Sgr A* under assumption of 
an isotropic stellar cusp). If the angular momentum transfer 
mediated by stars is very fast, then the surface column depth 
of the disk is much smaller (for the same M = 2-kR'Evr — 
const) than it is for the standard disk. If it becomes opti- 
cally thin, the disk will be much hotter for the same accre- 
tion rate. In fact, if it is unable to cool, the disk switches 
to a hot accretion mode, opening up an inner hole in the 
optically thick geometrically thin disk 2 . 

Due to interactions with the disk, embedded stars have 
Keplerian circular velocities. If two such stars meet in the 
disk, their relative velocity is small. In contrast stars out of 
the plane are expected to have a broad velocity distribution 
(e.g. Bahcall & Wolf 1976) and are prone to self-destruction 
in hard collisions (e.g. see §3 of Murthy et al. 1991, and 
references therein). Such a collision may involve two highly 
inclined stars or one star embedded in the disk and the other 
at a highly inclined orbit. If such an event happens in the 
inner disk region, and especially if the mass of the inner disk 
is not too much greater than the combined mass of the two 

2 Interestingly, this remark is also relevant to accretion disks 
around stellar-mass compact objects interacting with a cloud of 
comets - see Pineault & Landry 1994 and references therein 
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stars, then the resulting energy release will be sufficient to 
significantly speed up the accretion of the inner disk. The 
inner disk may then be quickly emptied out, and while it is 
being refilled the disk will radiate as a cut-off disk. 

As did the effects discussed in H4.4I these two mech- 
anisms for the inner disk removal are also of a temporary 
nature. Unless being flat as an accretion disk, the inner stel- 
lar cusp will self-destruct by binary collisions in a relatively 
short time, e.g. 10 7 — 10 s years for typical cluster param- 
eters (see, e.g., Fig. 6 of Murthy et al. 1991). Hence these 
mechanisms for the inner disk removal can be relevant to 
LLAGN only if these are relatively short lived phenomena 
in the SMBH/AGN activity cycles. 



6 DISCUSSION AND CONCLUSIONS 

Here we suggested that close stars (R <C 1 pc from the 
SMBH) created in the accretion disk by self-gravity in the 
bright quasar phase, may be present in the LLAGN disks. 
We found ( that it is unlikely that these disk stars signif- 
icantly contribute to the LLAGN optical/UV spectrum di- 
rectly. Hot young massive stars are ruled out because most 
LLAGN are weak in the UV; and low mass stars would have 
to be too numerous. Stellar contribution in the infra-red 
range is plausible, however, because SED of standard accre- 
tion disks are weak in that region. 

The effects of stars' mass, i.e. inertia, are much more 
realistic because they require a "small" number of low-mass 
stars embedded in the disk (say 10 2 — 10 4 stars). These stars 
can consume the disk faster than it can accrete on the SMBH 
(if no radial gap forms in the disk). In this case little gas may 
be flowing into the central SMBH. However, when fitting the 
optical/UV spectrum of a "typical" LLAGN, M81, we found 
that this model would only work if the disk were quickly 
and abruptly consumed by the embedded stars in the inner 
~ few hundreds Schwarzschild radii. Further work is needed 
to determine whether this is natural and may be sustained 
for long enough times to explain the LLAGN phenomenon. 

If a radial gap does form, then stars embedded in the 
disk do not accrete gas from the disk. However their presence 
will be still felt by the disk if the disk mass interior to the 
star is smaller than the mass of the star (or group of stars) . 
This appears to be especially important in the quiescent 
phase of limit-cycle unstable disks, in which case the stars 
will clog-up the inner accretion flow for time necessary for 
accumulation of enough disk matter to push the star in the 
disk. Since the disk interior to the star's location is quickly 
emptied out (Syer & Clarke 1995), such a star-terminated 
disk would appear as a cut-off disk. Finally, we noted that 
hard stellar collisions in the inner disk region may disrupt 
the disk - create temporary inner holes. 

We shall however notice that it is unlikely that the star- 
disk interactions discussed in this paper are the only expla- 
nation to the LLAGN "inner holes" in the accretion disk. As 
we have seen in H4.4I the character of star-disk interactions 
depends strongly on the mass of the SMBH. Observations 
span a very large range in the SMBH mass, i.e. almost 10 3 
(see e.g. Ho 1999, Ptak et al. 2004). Moreover, the accreting 
(stellar mass) binary systems may have inner holes in their 
accretion disks in the low luminosity states as well (e.g. see 
Esin, McClintock & Narayan 1997; McClintock et al. 2003). 



Therefore the star-disk interactions may be an additional 
reason for a complicated time-variability in a disk that is 
already unstable (e.g. Lin & Shields 1986; Siemiginowska et 
al. 1996) or time-dependent because of a changing supply 
of fuel at large distances, but they are unlikely to be the 
primary reason for such variations in the disk structure. 

Nevertheless with observations becoming progressively 
better, we should expect many surprises, especially exotic 
variability, from the effects of star-disk interactions. 
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